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Abstract 

QQ We analyze the lattice SU(2) Yang-Mills theory in Coulomb gauge. We show that the 

o 

^^ static gluon propagator is multiplicative renormalizable and takes the simple form D(|p|) = 

- X y^|p|2 + M^/\p\'^, proposed by Gribov through heuristic arguments many years ago. We find 

D ^ 

Q M = 0.88(l)GeV ~ 2y/a. 
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I. INTRODUCTION 

Non-Abelian gauge theories in Coulomb gauge have recently been the subject of 
extensive research. In particular, estimates for the ground state wave function have 
been derived through variational methods [H I2]- A key ingredient in such procedure 
is the static transverse gluon propagator D{\p\), which turns out to be compatible 
with the Gribov-Zwanziger scenario [31 E] . Unfortunately up to now ab initio lattice 
calculations have failed to deliver a conclusive comparable result. A first study carried 
out for SU(2) at fixed /3 = 2.2 [5l [6] indicated compatibility with Gribov's formula |3] 
in the IR but was inconclusive in the UV. Later studies in SU(2) and SU(3) showed 
strong scaling violations and contradictory IR results ior D^pl) [H |Hl E]- AH featured a 
UV behaviour at odds with simple dimensional arguments. Except for |S1IS], where the 
pole structure of the full propagator was analyzed, all these studies [3 El [9] have simply 
calculated the static propagator relying on the fact that once Coulomb gauge is fixed 
-^(IpI) = J2po^i\p\yPo) does not depend on the residual gauge freedom of temporal 
fields. One takes, however, de facto for granted that multiplicative renormalizability 
holds for the full propagator, although perturbative results seem to point at a more 
complex picture [TU], and that cutoff effects do not influence the result. To verify 
these assumptions we will fix analytically the residual temporal gauge and study the 
renormalization of the full propagator D{\p\,po). We will show that: i) the latter is 
not multiplicatively renormalizable; ii) the static propagator -D(|p|) is renormalizable 
only in the limit of continuous time, i.e. in the lattice Hamiltonian formulation. The 
resulting -D(|p|) = J_°° dpoD{\p\,pQ) satisfies the Gribov formula for all momenta [3]. 

II. SETUP 

We discretize the SU(2) Yang-Mills theory through the Wilson plaquette action 
S = l3J2xuiy(^ ~ ^ Tr f/^j,(x)), employing the standard Creutz algorithm (heath-bath 
plus over-relaxation). We take 2.15 < (3 < 2.6 and lattices L=12, 16, 20 and 24, with 
200 to 1000 measurements, depending on L. The lattice gluon propagator D'^Hp) is 



the Fourier transform of the gluon two-point function: 
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DSip) = {At{k)A%-k))^^ = 6^%,D{p) , p, = p{k,) = -sin ( V^ ) . (1) 



Here ^°(A;) is the Fourier transform of the lattice gauge potential A?^(x + /i/2) 
(see Eq. (|3])), a is the lattice spacing, p = {p,po) denotes the four-momentum and 
kfj, G (— L/2, +L/2] are the integer-valued lattice momenta. Results for the temporal 
component Dqq{p) will be presented in a forthcoming paper [llj. The Coulomb gauge 
is fixed for each configuration maximizing for every time slice 

^sW = ^E^r ^^(^'^) ' U!ix,t)=gix) U,{x,t) g^ix + l) (2) 

with respect to local gauge transformations g{x) G SU(2). The local maxima of Fg{t) 
satisfy for each t the differential lattice Coulomb gauge transversality condition for the 
gauge potentials: 

diAf{x,t)=Af{x + l/2,t)-Af{x-i/2,t) = 0, A^{x + fi/2) = ^(f/^W - f/^(x)t) . 

(3) 
In detail we: i) apply for the whole lattice a flip preconditioning [121 [13] as in Ref. [8]; 
ii) maximize Fg{t) via several standard over-relaxation and/or simulated annealing 
cycles after a random gauge copy for each time separately; iii) choose at every slice 
the copy with the highest Fg{t); iv) perform a random gauge transformation on the 
whole lattice and repeat from step i). The number of total recursions varies between 
10-40, depending on (3 and L. The copy of the whole lattice with the best functional 
value J2t ^gi^) ^^ chosen to select the gauge- flxed configuration. This procedure remains 
relatively efficient for every volume at higher (3, but in the strong coupling (low (3) region 
Gribov noise is no longer under control as the lattice size increases. This prevents us 
from going beyond L = 24. Future studies on larger lattices should rather adapt the 
improved algorithm proposed in Ref. [13] . 

Once the Coulomb gauge has been fixed, only space independent gauge transforma- 
tions g{t) can be performed at each time slice. In the continuum and on the lattice 



one can choose: 

= / d^xdf,A^{x,t) =^ do (fxAo{x,t) = <^ u(t) = —'^Uo{x,t) -^ const. 

X 

(4) 
Define 'u(t) = u(t)/Det [u(t)] G SU(2). Periodic boundary conditions make Tr n^-u(t) = 
TrP invariant under g{t). Eq. Q can thus be fixed recursively through g(t)u{t)g'^(t + 
I) = u = p^/L^ up to a g{0) satisfying [g{0),P] = 0. We can choose g{0) = 1. Of 
course u ^ 1 only for L -^ oo. We find 1 — Tr-u ~ 10"'^ already for L = 24, with 
leading corrections oc L^^. 

We always select spatial momenta indices k satisfying a cylinder cut P, [H] to 
minimize violations of rotational invariance. Time-like momenta are unconstrained. 
The scale is set re-expressing the lattice spacing a(/3) in units of the string tension 
a = (440MeV)2 [13]. 

III. RESULTS 

We find that the lattice bare propagator Df^dpljPo) satisfies over the whole p range: 

n f\^ \ fp(\p\) 9fi{z) po 

^/^(b1'Po) = ^^^^^ z=j^, (5) 

where |p|^(l + z^) is explicitly factorized since for large p"^ one expects D[^{\p\,po) ~ 
(IpI^+Po)^^ on dimensional grounds. fj3 and g/3 are dimensionless functions and without 
loss of generality we can choose (7/3(0) = 1. From Eq. ([s]) it is clear that a dynamical 
mass generation must occur in /^, i.e. it must be a function of \p\/M. We can check 
that taking J2po^p(\p\^Po) reproduces the results of Ref. |H]. The data for g/si^z) = 
{1 + z'^)Df){\p\,pQ)/D/j{\p\, 0) are shown in Fig. QJfor L = 24. Their leading behaviour is 
well described by a power law (1 + z^)". Sub-leading corrections can be parameterized 
equally well through different Ansatze, with all fits giving x^/d.o.f. between 0.8 and 
1.2. Without further theoretical input we have no reason to prefer one upon the other. 
For P > 2.3 the functions (7/3 vary quite consistently with L. In Table [l] we give 
our best estimates for a(/5) in the large L limit, assuming again a leading correction 
oc L^^. Within one or two standard deviations all values of a are compatible with 
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FIG. 1: Data for g^i^z) vs 1 + z^ in log-log scale, L = 24. For sake of readability not all /? are 
shown. 



1, i.e. D/3(|p|,po) might be po independent in the thermodynamic hniit. A study on 
larger, perhaps anisotropic lattices will be needed to clarify the matter. It is anyway 
clear from Eq. ([s]) that, as long as for our choices of /3 and L the functions gp{^z) 
remain different, as they do, D/3(|p|,po) (^o,'^ i^ot be multiplicative renormalizable. On 
the other hand, let us consider different lattice cut-offs for space and time ttg/at = 



/? 


2.15 


2.2 


2.25 


2.3 


2.4 


2.5 


2.6 


a(/3) 


1.01(2) 


1.02(3) 


1.02(3) 


0.99(5) 


0.96(6) 


0.8(1) 


0.6(2) 


Z{P,fj,) 


0.403(4) 


0.430(5) 


0.456(6) 


0.475(7) 


0.515(10) 


0.536(12) 


0.548(13) 



TABLE I: L ^ oo exponents a for the functions gpiz); renormalization coefficients Z(f3,^) 
as in Eq. dsl). 



^ > 1 [16]. For ^ — * oo (lattice Hamiltonian limit) D/^{\p\) oc |p1^^//3(|p|) up to the 
(possibly diverging) constant J dzgp{z)/{l + z^) and multiplicative renormalizability 
relies solely on //3(|p|). Eq. ^ also explains the scaling violations in Ref. [UlSlin]. Let 
us assume that //3(|p|) is renormalizable, define p = as\p\ and for simplicity neglect 
sub-leading corrections in gpiz). For large L we can approximate the discrete sum over 
Po by an integral and: 
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p 2tt 

where B{z,a,b) is the incomplete beta function. For $, ^ oo I explicitly goes to a 
\p\ independent constant, 5(1, 1/2, —a + 1/2), while in a standard lattice formulation, 
where ^ = 1, there is no way to avoid the extra |p| dependence I{2/p,a), which will 
violate multiplicative renormalizability. Indeed such function times the tree level result 
\p\^^ compares well with the behaviour observed in [3 [HI |9]. 

From the above discussion it is clear that the po dependence is immaterial to the 
static propagator in the Hamiltonian limit, all the relevant information being encoded 
in fi3{\p\)- Since the limit ,^ — > oo is out of the reach of available simulations, to extract 
fp from simulations at standard sizes and coupling (i.e. spacing) one could use only 
results at fixed po, e.g. //3(|p|) = -D/3(|p|,0), at the cost of loosing data. Although this 
works, we propose here an alternative procedure which makes use of the results for all 
energies, increasing the statistics. We divide the z dependence out of the data using 



the fitted functional form for (7/3(2;) and average the result over po: 

The factor |p|~^ only serves to restores dimensions when comparing with the Hamilto- 
nian result, according to Eq. ([6]), and is immaterial to renormalization. To minimize 
finite size effects we exclude spatial momenta not satisfying a cone cut [21 HI] for all 
lattices with (a{P)L)^ < (2.5fm)^. The Dp{\p\) calculated form Eq. Q are multiplica- 
tive renormalizable via D^(|p|)| = Z{(3, ^)Dj3{\p\), where fi defines the renormalization 
point. We choose it such that (see Table [l]): 

D,m\,=oo = Zif3,^^)\,=^Dpm\ ,r- 7^, (8) 

|p|->oo Z\p\ 

corresponding to the naive UV continuum behaviour. Fig. [2] shows the result of our 
procedure. We have performed fits of D^(|p|) allowing at the same time a power 
law in the IR, \'p\'^/M'^~'^, and a power law plus logarithmic corrections in the UV, 
Ipl-'^lloglpl"'. We get q = 0.99(1), r = 1.002(3), s = 0.002(2) and M = 0.88(1) GeV, 
with x^/d.o.f. all in the range 3.2-3.3, in agreement with the UV and IR analysis in 
Ref. [21 [T71 [18]. Since our data show excellent IR-UV symmetry for p ^->- M'^/p (see 
Fig. [I]), we constrain g = r = l,s = and fit the whole result through the Gribov 
formula 

DM)\,=oo = I ^ „^, (9) 

We find just as good agreement (x^/d.o.f. = 3.3) again with M = 0.88(l)GeV ~ 2^/(x. 
This compares well with M ^ 0.5 — 0.8 GeV from the IR analysis in Ref. [51 Ej, 
given that only j3 = 2.2 was used there. Apart from unavoidable residual violations of 
rotational invariance, the highest noise in the data occurs around the peak, \p\ ~ M. 
Only better gauge fixing and larger volumes can help here. We have also collected data 
for 2.6 < (3 < 2.8. They are still in perfect agreement with the above formula, however 
since physical volumes are very small we cannot claim to have control over systematic 
effects, so we omit them in the above fits. An investigation on higher volumes would 
be welcome. 
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FIG. 2: Gluon propagator. Circles are (a few) data for \p\ — > M^/|p|, the continuum line is 
the fit to Gribov's formula and the dot-dashed line the result of the Hamiltonian approach 
of Ref. [H]. 

IV. CONCLUSIONS 

We have shown that on the lattice the static gluon propagator in Coulomb gauge is 
multiplicative renormalizable only in the Hamiltonian limit. This calls for analogous 



investigations in the continuum, where to our knowledge renormahzabihty has not 
been proven in Coulomb gauge. We have also given a procedure to extract the relevant 
information from the data at accessible lattice sizes and spacings. Summarizing this 
procedure: i) calculate the full bare propagator D^(|p|,po) after fixing the residual 
gauge as in Eq.Q; ii) fit to the data the function gp{z) arising from the factorization in 
Eq.([5]); iii) divide out g^iz) from D^(|p|,po) to isolate itspo independent part ff3{\p\); iv) 
define D/3(|p|) = //3(|p|)/|p|. The so obtained static propagator Dp^pl) is multiplicative 
renormalizable, UV-IR symmetric and is well described by Gribov's formula over the 
whole momentum range. Its infrared and ultraviolet behaviours are in good agreement 
with the results obtained in the variational approach to continuum Yang-Mills theory 
in Coulomb gauge ^ . Both the infrared analysis of the Dyson-Schwinger equations in 
this approach and its numerical solution yield for the infrared exponent g = 1 [T71IT8] . 
We will compare further quantities obtained in the variational approach [2] with the 
lattice using the procedure proposed above [TT] . 
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